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Abstract. Models of thermal emission of neutron stars, 
presumably formed in their atmospheres, are needed to 
infer the surface temperatures, magnetic fields, chemical 
composition, and neutron star masses and radii from the 
observational data. This information, supplemented with 
model equations of state and neutron star cooling mod- 
els, is expected to move us further in understanding the 
fundamental properties of the superdense matter in the 
neutron star interiors. The neutron star atmospheres are 
very different from those of usual stars due to the immense 
gravity and huge magnetic fields. In this presentation we 
review the current status of the neutron star atmosphere 
modeling and present most important results. 
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1. Introduction 

A systematic study of X-ray emission from isolated neu- 
tron stars (NSs), including radio pulsars, has started af- 
ter the launch of the Einstein and EXOSAT space obser- 
vatories. These studies have shown that, generally, there 
are two different components of the NS X-ray emission 
- thermal and nonthermal. The nonthermal component 
with a power-law spectrum, observed from many radio 
pulsars, is believed to originate from the pulsar's magne- 
tosphere, while the thermal component is emitted from the 
NS surface layers (atmospheres) . The thermal radiation is 
particularly interesting because it can provide important 
information about the NS: its surface temperature, mag- 
netic field, and chemical composition, as well as the NS 
radius and mass. Measuring these parameters for a sample 
of NSs is necessary for studying the thermal evolution of 
NSs and constraining the equation of state and composi- 
tion of the superdense matter in the NS interiors (see the 
review by Yakovlev et al. in these Proceedings). 

Thermal X-ray emission from the NS surface had been 
discussed by Chiu & Salpeter and Tsuruta in 1964, before 
NSs were discovered, and well before their thermal emis- 
sion was detected with Einstein and EXOSAT (e.g., Cheng 
& Helfand 1983; Brinkmann & Ogelman 1987; Cordova et 
al. 1989; Kellet at al. 1987). Many new observational re- 
sults on the NS thermal radiation were obtained in 1990' 



with the ROSAT, ASCA, and EUVE satellites (see Becker 
& Pavlov 2002 for a review) . Thermal radiation from a few 
NSs was also detected in the optical-UV energy range with 
the Hubble Space Telescope (e. g., Pavlov, Stringfellow & 
Cordova 1996; Walter & Matthews 1997). Currently op- 
erating X-ray observatories, Chandra and XMM-Newton, 
are providing new excellent data on X-rays from NSs (see 
the contributions by Weisskopf, Becker, and Pavlov, Za- 
vlin & Sanwal in these Proceedings). 

To interpret these observations, one needs reliable 
models for the NS thermal radiation. The fact that the 
spectrum of radiation emergent from a NS atmosphere 
can be very different from a blackbody spectrum, and 
the angular distribution can be very far from isotropic, 
particular in a strong magnetic field, has been recognized 
long ago. For instance, Pavlov & Shibanov (1978) calcu- 
lated spectra and angular distributions of radiation from a 
strongly magnetized NS atmosphere assuming the source 
function grows inward linearly in the emitting layers. First 
self-consistent models for NS atmospheres were developed 
by Romani (1987), for low magnetic fields, and Shibanov 
et al. (1992), for strong magnetic fields. Since then, var- 
ious aspects of the NS atmosphere modeling have been 
investigated in many papers. Below we will overview the 
current status of this field and summarize some important 
results. 



2. Distinctive features of NS atmospheres 

The NS thermal radiation emerges, as in usual stars, in 
superficial layers which can be in a gaseous state (atmo- 
sphere) or a condensed state (liquid or solid surface), de- 
pending on surface temperature, magnetic field and chem- 
ical composition. A condensed surface forms at low tem- 
peratures and very strong magnetic fields. For instance, 
according to the estimates by Lai & Salpeter (1997), hy- 
drogen is condensed in surface layers if T & 1 x 10 5 K 
at B = 1 x 10 13 G (T < 1 x 10 6 K at B = 5 x 10 14 
G). At higher temperatures and/or lower magnetic fields, 
hydrogen does not condensate and forms an atmosphere. 

In both usual stars and NSs the properties of radi- 
ation emergent from an atmosphere strongly depend on 
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its chemical composition. The difference between usual 
stars and NSs is that in NSs we expect the emitting layers 
to be comprised of just one, lightest available, chemical 
element, not a mixture of elements, because heavier el- 
ements sink into deeper layers due to the immense NS 
gravitation (Alcock & Illarionov 1980). For instance, even 
a small amount of hydrogen (with a surface density of 
~ 1 g cm" 2 ) is sufficient for the radiation to be indis- 
tinguishable from that emitted from a purely hydrogen 
atmosphere. Such an amount of hydrogen, > 1O _2O M0, 
can be delivered onto the NS surface by, e. g., accretion 
from the interstellar medium and/or fallback of a fraction 
of the envelope ejected during the supernova explosion. 
If no hydrogen is present at the surface, a heavier chem- 
ical element is responsible for the radiative properties of 
the NS atmosphere. However, a mixture of elements can 
be observed in the emitting layers if the NS is experienc- 
ing accretion with such a rate that the accreting matter 
is supplied faster than the gravitational separation occurs. 
Therefore, it is important to construct the NS atmosphere 
models for a variety of surface compositions, both 'pure' 
and mixed. 

The gravitational separation of elements is one of the 
consequences of the enormous gravity at the NS surface, 
with the gravitational acceleration g ~ 10 14 — 10 15 cm s~ 2 . 
The gravity makes NS atmospheres very thin, ~ 0.1 — 10 
cm, and dense, p ~ 10~ 2 — 10 2 g cm" 3 . Such a high den- 
sity causes strong nonideality effects (pressure ionization, 
smoothed spectral features) which must be taken into ac- 
count (e. g., Pavlov et al. 1995). In addition, the strong 
gravitational field bends the photon trajectories near the 
NS surface (Pechenick, Ftaclas & Cohen 1983), as illus- 
trated by the sketch in Figure |l|. This effect depends on 
the gravitational parameter g T = (1 — 2GM/c 2 R) 1 ^ 2 , and 
it can even make the whole NS surface visible if the NS 
is massive enough, M > 1.8 (lOkm/i?) M Q . In particu- 
lar, the gravitational bending strongly affects the observed 
pulsations of thermal emission (Zavlin, Shibanov & Pavlov 
1995a). 

Huge magnetic fields, B ~ 10 11 — 10 14 G, in the sur- 
face layers of many NSs change the properties of the at- 
mospheric matter and the emergent radiation even more 
drastically. Strongly magnetized NS atmospheres are es- 
sentially anisotropic, with radiative opacities depending 
on the magnetic field and the direction and polarization 
of radiation. Moreover, since the ratio of the cyclotron 
energy, E cc = HeB/m c c, to the Coulomb energy is very 
large (e.g., 7 = E cc /(1 Ry) = 850 {B/10 12 G) for a hydro- 
gen atom), the structure of atoms is strongly distorted by 
the magnetic field. For instance, the binding (ionization) 
energies of atoms are increased by a factor of ~ In 2 7 (e. g., 
the ionization potential of a hydrogen atom is about 310 
eV at B = 10 13 G). This, in turn, significantly modifies 
ionization equilibrium of the NS atmospheric matter. An- 
other important effect is that the heat conductivity of the 
NS crust is anisotropic (it is higher along the magnetic 




Fig. 1. Bending of photon trajectories in a strong gravi- 
tational field near the surface of a NS with rotational and 
magnetic axes ft and m, respectively. 

field). This results in a nonuniform surface temperature 
distribution (Greenstein & Hartke 1983), which leads to 
pulsations of the thermal radiation due to NS rotation. 

3. Atmosphere models with low magnetic fields 

Some NSs (e.g., millisecond pulsars) have relatively low 
magnetic fields, B 10 9 G. In such magnetic fields the 
electron cyclotron energy, E cc Ss 0.01 keV, is lower than 
the binding energy of atoms and thermal energy of par- 
ticles. As a result, the effect of the field on the radiative 
opacities and emitted spectra is negligible, at X-ray en- 
ergies, E il 0.1 keV. Hence, to model low-field NS at- 
mospheres, one can merely put B = 0. These models are 
applicable to millisecond pulsars, NS transients in quies- 
cence (e.g., Rutledge et al. 2002), and, perhaps, to some 
radio-quiet isolated NSs. One should remember, however, 
that the magnetic field effects on spectral opacity at opti- 
cal wavelengths must be included even for these low fields 
if the models are applied for interpreting the optical spec- 
tra. 

First models of low-field NS atmospheres were pre- 
sented in the pioneering work by Romani (1987). Since 
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Fig. 2. Radiative spectral opacity fcj, for hydrogen, helium 
and iron plasmas at different temperatures and densities 
(magnetic held B = 0). The iron opacity is provided by 
the OPAL group (Rogers & Iglesias 1994). 



then, models for various surface chemical compositions 
have been developed by Rajagopal & Romani (1996), Za- 
vlin, Pavlov & Shibanov (1996), Pavlov & Zavlin (2000a), 
Werner & Deetjen (2000), Pons et al. (2002), Gansicke, 
Braje & Romani (2002). 

3.1. General approach 

As the thickness of a NS atmosphere is much smaller than 
the NS radius, R f» 10 km, NS atmospheres can be con- 
sidered in the plane-parallel approximation. In addition, 
rather high densities of the surface layers allow one to 
consider the atmospheres to be in local thermodynamic 
equilibrium. 

The standard approach for the atmosphere modeling 
includes solving of three main equations. The hrst one 
is the radiative transfer equation for the specific spectral 
intensity I v (e. g., Mihalas 1978): 
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Fig. 3. Temperature dependences T(y) in nonmagnetic 
atmosphere models with different effective temperatures 
and chemical compositions. 

where fi is cosine of the angle between the normal to 
the atmosphere and the wave-vector of radiation, y is the 
column density (dy = pdz, with z being the geometri- 
cal depth), fc„ = a v + a v is the total radiative opacity 
which includes the 'true absorption' (a„) and scattering 
(<r„) opacities, S v = (a v J v -I- a v B v )k^ is the source func- 
tion, Jv = h J, Ivdfi is the mean spectral intensity, and 
B v is the Planck function. The boundary condition for 
this equation is I v = for fi < at y — 0, assuming no 
incident radiation at the NS surface. 

NS atmospheres are usually considered to be in ra- 
diative and hydrostatic equilibrium. The hrst condition 
implies that the total energy flux through the atmosphere 
is constant and transferred solely by radiation, 



fj,I v dfi 



(2) 



where T e g is the effective temperature, and <tsb is the 
Stefan-Boltzmann constant. The second condition means 
that the atmospheric pressure is P = g y (the radiative 
force is insignificant unless T e g <; 10 7 K). Finally, these 
three equations are supplemented with the equation of 
state for the atmospheric plasma and equations of ioniza- 
tion equilibrium. The latter is needed for computing the 
electron number density and the fractions of ions in dif- 
ferent stages of ionization to obtain the radiative opacity 
with account for free-free, bound-free and bound-bound 
transitions (see examples in Fig. |^)Q A detailed descrip- 

1 Opacities of heavy chemical elements used in this work 

were calculated by the OPAL group: 

|http : //www-phys . llnl . gov/Research/OPAL/ index .html 
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Fig. 4. Spectral fluxes of emergent radiation in nonmagnetic atmosphere models for several effective temperatures and 
chemical compositions. 



tion of the model computations can be found, e.g., in Za- 
vlin ct al. (1996). 



3.2. Results 

To explain the main properties of the atmosphere emis- 
sion, we first present in Figure ^ the most important prop- 
erty of the atmospheric structure responsible for the emer- 
gent radiation — the depth dependence of the tempera- 
ture in the NS surface layers, T(y). 

Figure f| presents the spectral fluxes of emergent radia- 
tion at a local surface point, F u — J Q /i/„d// (at y = 0), for 
several effective temperatures and chemical compositions 
(pure hydrogen, helium and iron, and a solar mixture), to- 
gether with the blackbody spectral fluxes at the same T e g. 
The model spectra differ substantially from the blackbody 
spectra, particularly in high-energy tails of the spectra of 
light-element (hydrogen and helium) atmospheres^. The 
reason for such behavior is in the rapid decrease of the 
light-element opacities with energy (~ E~ 3 , see the upper 
and middle panels in Fig. 0), so that the high-energy radi- 



2 Th e nonmagnetic hydrogen atmosphere models are avai l- 
able at tittp : //legacy .gsfc .nasa.gov/docs/xanadu/xspec 



ation is formed in deeper and hotter layers (Fig. |3|). Since 
the radiative opacity is large at lower energies, the low- 
energy photons escape from very superficial layers with a 
temperature T(0) < T e g. Therefore, the low-energy tails 
of the model spectra are suppressed with respect to the 
Ray leigh- Jeans spectra at the same T e ff. 

The spectra emitted from the heavy-element atmo- 
spheres (iron and solar mixture) exhibit numerous spectral 
lines and photoionization edges (e. g., M, L and K spec- 
tral complexes in the iron spectra) produced by ions in 
various ionization stages. Generally, they are closer to the 
blackbody radiation because the energy dependence of the 
heavy-element opacities is, on average, flatter than that 
for the light elements (see examples in the bottom panel 
of Fig. ||). However, local deviations from the blackbody 
spectra can be very substantial. Moreover, the spectral 
features can be rather strong even if metals contribute as 
little as 0.05% of the total NS atmosphere composition 
(see Fig. ||). 

Although the opacity of the atmospheric plasma is 
isotropic in the nonmagnetic case, the emitted radiation 
show substantial anisotropy, i. e., the specific intensity I v 
depends on the direction of emission due to the limb- 
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Fig. 5. Spectral fluxes of emergent radiation in nonmag- 
netic models with different metal abundances Z (Z = 
2.0% corresponds to the standard solar mixture of ele- 
ments — see Grevesse & Noels 1993). 



darkening effect (Fig. At large values of the angle 9 
between the normal to the surface and the wave-vector 
the emerging photons are produced in shallow layers with 
lower temperatures. The anisotropy depends on photon 
energy and chemical composition of the atmosphere. This 
effect should be taken into account to model thermal ra- 
diation from a nonuniform NS surface (see Section 5). 

The emergent radiation depends also on the surface 
gravity: a stronger gravitational acceleration increases the 
density of the atmospheric plasma and enhances the non- 
ideality effects, which results in weaker (more smoothed) 
spectral features (see examples in Fig. 6 of Zavlin et al. 
1996). However, this effect is rather small and might be im- 
portant only for analyzing observational data of extremely 
good statistics. 

4. Atmosphere models with strong magnetic fields 

It is commonly believed that most NSs (at least ac- 
tive radio/X-ray/7-ray pulsars) have rather strong sur- 
face magnetic fields, B ~ 10 10 -10 14 G. Such fields dras- 
tically change the properties of NS atmospheres, as dis- 
cussed above. Magnetic hydrogen models have been de- 
veloped by Shibanov et al. (1992), Pavlov et al. (1994), 
Shibanov & Zavlin (1995), Pavlov et al. (1995), and Za- 
vlin et al. (1995b). These models use simplified radiative 
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Fig. 6. Spectra of specific intensities l v emitted from an 
iron atmosphere with T e g = 1 x 10 6 K, at different angles 
9 (/* = cos 9). 



opacities of strongly magnetized, partially ionized plasma, 
which do not include the bound-bound transitions. How- 
ever, they are considered to be reliable enough in the case 
of high temperatures (T e g ^ 10 6 K at typical radio-pulsar 
fields ~ 10 12 G), when the atmospheric plasma is almost 
fully ionized even in the strong magnetic fields. Recently, 
completely ionized hydrogen models for superstrong mag- 
netic fields, B ~ 10 14 ~10 15 G, have been presented in 
a number of papers (Bezchastnov et al. 1996; Bulik & 
Miller 1997; Ozel 2001; Ho & Lai 2001; Zane et al. 2001), 
concerned mainly with the vacuum polarization effects 
(Pavlov & Gnedin 1984) and the ion cyclotron lines whose 
energies get into the X-ray band at B <; 2 x 10 13 G. A set 
of magnetic iron models was constructed by Rajagopal, 
Romani & Miller (1997), with the use of rather crude ap- 
proximations for the very complicated properties of iron 
ions in strong magnetic fields. The iron atmosphere spec- 
tra show many prominent spectral features which could 
be very useful to measure the NS magnetic field and the 
mass-to-radius ratio. 

General approach 

Modeling magnetic NS atmospheres is generally similar 
to the nonmagnetic case. The main difference is that the 
atmospheric radiation is polarized, and the radiative opac- 
ities depend on the polarization and direction of radiation. 
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Fig. 7. Spectral fluxes of radiation emerging from magnetic hydrogen atmospheres, for several effective temperatures 
and magnetic fields perpendicular to the surface. 



Gnedin & Pavlov (1974) described the radiative transfer 
in a strongly magnetized plasma in terms of coupled equa- 
tions for specific intensities of two normal modes, I v _\ and 
Iv,2, with different polarizations and opacities: 

fji—Iv,j(n) = k Utj (n)I Utj (n)- 



2 

E 

,»=i 



j> dn' I v ,i(n') a U}ij (n' ' ,n) 



(3) 



where n is the (unit) wave-vector, a„j is the absorption 
opacity for the j-th mode, and a u ,ij is the scattering opac- 
ity from mode i to mode j. It should be noted that the 
opacity depends on the angle between the wave- vector and 
the magnetic field, so that /„ depends not only on 9, but 
also on 9b, the angle between the local magnetic field and 
the normal to the NS surface. Similar to the nonmagnetic 
case, equations (3) are supplemented with the equations 
of radiative and hydrostatic equilibrium. To overcome the 
problems with the sharp angular dependence of the ra- 
diative opacities (Kaminker, Pavlov & Shibanov 1982), 
a two-step method for modeling of magnetic NS atmo- 
spheres was developed (Pavlov et al. 1994; Shibanov & 



Zavlin 1995). At the first step, the radiative transfer is 
solved in the diffusion approximation for the mean inten- 
sities J u j = (Ait)^ 1 § I V j(n) dn: 



_d, ± T _- 

dy v ' j dy v ' j ~ " 



2 



(4) 



where d U j is the diffusion coefficient, 



a v j and cf„ are the 
angle-averaged absorption and scattering opacities (see 
Pavlov et al. 1995 for details). Next, the atmospheric 
structure obtained at the first step is corrected using an 
iterative procedure applied to the exact equations of the 
radiative transfer. 

4-2. Results 

Figure [j] shows polarization-summed spectral fluxes of the 
emergent radiation, F v = J Q jU(J M) i +7 /tl 2)d/x, emitted by 
a local element of the NS surface, for several values of 
effective temperature and magnetic field perpendicular to 
the surface (8b = 0). The main result is that the magnetic 
model spectra are harder than the blackbody radiation of 
the same T e g, although they are softer than the low-field 
spectra. Similar to the low-field case, this is explained by 
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Fig. 8. Spectral flux of radiation emerging from an atmo- 
sphere with B = 10 14 G, with a strong proton cyclotron 
line. 

the temperature growth with depth and the opacity de- 
crease at higher energies, which is more gradual (oc E^ 1 
for the mode with smaller opacity) compared to the non- 
magnetic case. At lower effective temperatures, T e ff 10 6 
K, the photoionization opacity is important, and the de- 
pendence of the spectra on magnetic field becomes more 
pronounced (see the right panel in Fig. ^). The field de- 
pendence is most clearly seen, at T e g sufficiently high to 
provide strong ionization, in the proton cyclotron lines 
centered at energies E cp = 6.3 (B/10 12 G) eV. If the mag- 
netic field is very large, B > 10 14 G, the proton cyclotron 
line shifts into the X-ray band and can be very strong (see 
Fig. ||). On the other hand, if the magnetic field is not so 
large, B — 10 10 -10 12 G, the NS atmosphere spectra may 
exhibit the electron cyclotron lines in the X-ray band, at 
E cc = 11.6 (B/10 12 G) keV (see Fig. |). 

First calculations of hydrogen atmosphere models 
which include bound-bound transitions show that spec- 
tral lines, considerably broadened by the motional Stark 
effect (Pavlov & Meszaros 1993; Pavlov & Potekhin 1995) 
become prominent at T e g ^ 5 x 10 5 K. The strongest 
line corresponds to the transition between the ground 
state and the lowest excited state; its energy is E » 
[75(l + 0.131n(5/10 13 G) + 63(5/10 13 G)] eV, at B ~ 10 13 
G (see examples in Fig. [l(]). 

Radiation emerging from a magnetized NS atmosphere 
is strongly anisotropic. Angular dependences of the lo- 
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Fig. 9. Spectral flux of radiation emerging from an at- 
mosphere with B — 3 x 10 11 G, with a strong electron 
cyclotron line. 
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Fig. 10. Spectral fluxes of radiation emerging from hy- 
drogen atmospheres with B — 2.4 x 10 12 G and different 
values of T e g, with bound-bound transitions included. 
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Fig. 11. Angular dependences of polarization-summed 
specific intensities at different photon energies E in the 
hydrogen atmosphere model with T e g = 1 x fO 6 K and 
magnetic field B = 4 x f 12 G, perpendicular to the sur- 
face. 



cal specific intensities, l v = J„.i + I v ,2, show a compli- 
cated "pencil-plus-fan" structure — a narrow peak along 
the direction of the magnetic field (where the atmo- 
spheric plasma is most transparent for the radiation) , and 
a broader peak at intermediate angles. The widths and 
strengths of the peaks depend on magnetic field and pho- 
ton energy (see examples in Fig. [ll]). Obviously, it is very 
important to account for this anisotropy while modeling 
the radiation from a NS with nonuniform surface magnetic 
field, effective temperature, and chemical composition. 



5. Thermal radiation as seen by a distant 
observer 

Results presented in Sections 3 and 4 describe spectral ra- 
diation emitted by a local element at the NS surface. The 
effective temperature and/or magnetic field distributions 
over the NS surface can be nonuniform (for example, if 
a NS has a dipole magnetic filed, the effective tempera- 
ture decreases from the magnetic poles to the equator). 
To calculate the total NS emission, one has to integrate 
the local intensities, computed for local temperatures and 
magnetic fields, over the visible part of the NS surface S, 
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Fig. 12. Spectral fluxes of radiation from the whole sur- 
face of a nonmagnetic NS covered with an iron atmosphere 
at T c ff = 1 x 10 6 K. The (unredshifted) spectra are calcu- 
lated for different NS rotation periods P, at the inclination 
of rotational axis £ = 90°. 

with account for the gravitational redshift and bending of 
photon trajectories: 

F ( E ) = 5r-H2 / »I(E/9r)dS, (5) 
« s 

where D is the distance from the NS to the observer, 
and E is the observed (redshifted) energy^]. More details 
about the integration over the NS surface can be found in 
Pavlov & Zavlin (2000b). In should be noted that if a NS 
has a nonuniform (e.g., dipole) distribution of the mag- 
netic field, the integration broadens the spectral features 
(in particular, cyclotron resonance lines — see Zavlin et 
al. 1995b; Pavlov & Zavlin 2000b; Zane et al. 2001). 

If the NS is a fast rotator, one should take into account 
the Doppler shifts of energies of photons emitted from 
surface elements moving with different radial velocities. 
Maximum values of these velocities, v r = 2ttRP~ 1 sin£ 
(P is the rotation period, ( the inclination of the rotation 
axis), can be as high as 10%-15% of the speed of light 
for millisecond periods. Figure |l2| presents (unredshifted) 
spectra emitted from the whole NS surface covered with 

3 To take into account the interstellar absorption, a factor 
of exp[— riHcr c ff(£0] should be added in eq. (5) [a c fi(E) is the 
absorption cross section per hydrogen atom]. 
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Fig. 13. Energy-integrated light curves of a rotating NS with different distributions of surface temperature and 
magnetic field (from Shibanov ct al. 1995). The temperature is either uniform, T — T p , or nonuniform, T = T(#), 
as given by computations with account for the anisotropy of the heat conductivity in the NS crust ($ is magnetic 
colatitude, T p and B p are the effective temperature and the field at the magnetic poles). The NS rotational axis is 
perpendicular to both the magnetic axis and the line of sight (a — ( = 90°). 



an iron atmosphere (nonmagnetic case). For slowly rotat- 
ing NSs, P <; 20 ms, there is no significant effect on the 
emergent radiation. Fast rotation, P Ss 10 ms, may lead 
to complete smearing of weak and narrow spectral lines, 
provided sin£ is large enough, leaving only most promi- 
nent spectral jumps around the strongest photoionization 
edges. 

In the case of radio pulsars, a fraction of the observed 
thermal X-ray radiation is emitted from small hot spots 
(polar caps) around the NS magnetic poles. These polar 
caps can be heated up to X-ray temperatures by relativis- 
tic particles impinging on the poles from the acceleration 
zone in the pulsars' magnetosphere. In the case of mil- 
lisecond pulsars, with characteristic ages of 10 8 -10 9 yr, 
the whole NS surface (except for the polar caps) is be- 
lieved to be too cold (T c s 10 5 K) to be observable in 
X-rays, so the detected thermal emission can be modeled 
as emitted from small heated spots: 

F{E)=g t ^I{E/g T ,e*), (6) 

where the apparent spot area S a and the angle 9* between 
the wave-vector and the radius-vector to the hot spot 
are computed with account for the effect of gravitational 
bending. These quantities depend on the angles a (be- 
tween the rotational and magnetic axes) and £ (between 
the rotational axis and the line of sight — see Fig. |l|), and 
the gravitational parameter g r (see Zavlin et al. 1995a for 
details). 

The flux given by equations (5) or (6) varies with the 
period of NS rotation. One can obtain a large variety of 
pulse profiles at different assumptions on the angles a and 



£ and the NS mass-to-radius ratio. Examples of pulse pro- 
files computed for radiation from the whole NS surface 
are shown in Figure [l3|; pulse profiles of thermal radia- 
tion from heated polar caps are presented by Zavlin et al. 
(1995a). 

Although the model atmosphere spectra are different 
from the blackbody spectra, very often an observed ther- 
mal spectrum can be fitted equally well with the black- 
body and NS atmosphere models, particularly when the 
energy resolution is low and/or the energy band is nar- 
row. However, the parameters obtained from such fits are 
quite different, especially when the hydrogen or helium 
atmospheres are used. Since the light-element atmosphere 
spectra are much harder than the blackbody spectra at 
the same effective temperature, atmosphere model fits give 
temperatures T atm significantly lower than the blackbody 
temperature Tbt>, with a typical ratio Tbb/T atm ~ 2-3 (see, 
for example, Pavlov et al. 1996; Pons et al. 2002). On 
the other hand, to provide the same total energy flux, the 
blackbody fit gives a smaller normalization factor, propor- 
tional to S/D 2 , than the atmosphere model fit. In other 
words, the light-element atmosphere fit gives a consider- 
ably larger size of the emitting region, S^tm/Sbb ~ 50-200, 
for the same distance to the source. 

6. Conclusions 

In the recent decade, substantial progress has been made 
in modeling atmospheres of isolated NSs. Best inves- 
tigated cases are nonmagnetic atmospheres and fully- 
ionized light-element atmospheres with strong magnetic 
fields. The atmosphere models have been applied to the 
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interpretation of thermal emission from NSs of different 
types. For instance, Pavlov & Zavlin (1997), Zavlin & 
Pavlov (1998), and Zavlin ct al. (2002) analyzed the X-ray 
emission from the millisecond pulsar J0437 -4715 using the 
nonmagnetic light-element atmosphere models; Pavlov et 
al. (2001) applied the magnetic hydrogen models to the 
analysis of radiation from the Vela pulsar; Zavlin, Pavlov 
& Triimper (1998) and Zavlin, Triimper & Pavlov (1999) 
used the magnetic light-element models to interpret the 
thermal emission from NSs in the supernova remnants 
PKS 1209-51/52 and Puppis A. More examples are pre- 
sented in the contribution by Pavlov et al. in this volume. 

However, a number of problems in the atmosphere 
modeling remains to be solved. First of all, investigations 
of the structure of various atoms, molecules, and molec- 
ular chains in strong magnetic fields, as well as radiative 
transitions in these species (Pavlov 1998), are necessary to 
construct magnetic atmosphere models of different chem- 
ical compositions. First efforts in this direction are being 
undertaken (Mori & Hailey 2002). Particularly interesting 
are the (virtually unknown) radiative properties of mat- 
ter in superstrong magnetic fields, B il 10 14 G, apparently 
found in anomalous X-ray pulsars and soft gamma-ray re- 
peaters. Further work is needed on radiative properties of 
nonidcal plasmas and condensed matter. These investiga- 
tions will eventually result in more advanced models for 
thermal radiation of isolated NS, needed for the interpre- 
tation of high-quality observations of these objects with 
the Chandra and XMM-Newton X-ray observatories. 
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